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The superconductivity in cold dense QCD with three massless flavors is analyzed in the framework 
of the Schwinger-Dyson equation. The set of two coupled gap equations for the color antitriplet, 
flavor antitriplet (3,3) and the color sextet, flavor sextet (6,6) order parameters is derived. It is 
shown that the antitriplet-antitriplet gives the dominant contribution to the color-flavor locked order 
parameter, while the sextet-sextet is small but nonzero. 
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The idea of color superconductivity in dense QCD is not new B,0|. However, it is only now that we start to 
appreciate the potential relevance of this phenomenon for heavy ion collisions and for physics of neutron (or quark) 
stars. The recent progress was stimulated by the observation that the cold dense quark matter could undergo a phase 
transition into a color superconducting state with the value of the order parameter as large as 100 MeV ||||]. For 
comparison, the most optimistic estimates of Ref. were at the order of 1 MeV or less. Not surprisingly, the new 
estimate triggered the burst of studies on the subject mlq|. 

The approach of Refs. yjj] was based on the analysis of QCD in the approximation in which the interaction between 

quarks is modeled by instantons. Such an approach is best suited for studies of the low density quark matter. At 

^ ' densities relevant for the phenomenology of heavy ion collisions and physics of neutron stars, however, the instanton 

l/*~) . approach may just start to fail. Therefore, it was of equally great importance to rederive the same result in the 

framework of the weakly coupled QCD at asymptotically large densities. Along this direction, first it was pointed 

by Pisarski and Rischke and independently by Son H] that the unscreened gluon modes of magnetic type could 

provide the interaction necessary to generate a large value of the order parameter. Later, using the approach of the 

qs \ Schwinger-Dyson (SD) equation, it was shown that this hypothesis was indeed correct p|-^2|. 

Q\ • In this paper, we generalize the analysis of Refs. |9|-[l4| to the case of dense QCD with three massless flavors. In 

"^ \ particular, we pay special attention to the color-flavor structure of the order parameter. As is known, QCD with three 

massless flavors should reveal the color-flavor locked phase where the original SU(3)c x SU(3)r X SU(3)l X U{\)b 

symmetry of the action breaks down to SU(2>)c+l+r x Zi (iJJMI- From the symmetry arguments alone ||, the 

Qh, order parameter should consist of the color antitriplet, flavor antitriplet (3, 3) and the color sextet, flavor sextet (6, 6) 

contributions. This is also confirmed by the explicit solution in the effective theory obtained from the instanton liquid 

model. There, it is shown that, while the dominant contribution to the order parameter is the antitriplet-antitriplet 

k*" ■ one, the sextet-sextet admixture is still nonzero, although small |L4|. Below, we will demonstrate that the same 

conclusion remains true when the interaction is mediated by the perturbative one-gluon exchange. 

As in Ref. M, instead of working with Dirac spinors, it is more convenient to introduce the eight component 



where ip c — C-tp T and C is a unitary charge conjugation matrix, defined by C~ 1 7 /J C = — 7J and C = —C T . In this 
new notation, the bare fermion propagator is defined by 



(G(°)(p)) 



l = -i(* + ^° , ° „]. (2) 



i> - hj° 



where j> = 7"p„ and fi is the chemical potential. Regarding the notation, we also note that the fermion field and 
its propagator bear the color (a, b, . . .) and flavor (i, j, . . .) indices. Often, however, when it does not cause any 
confusion, we will omit them. 

The SD equation in the improved rainbow approximation reads H 
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Here / D llu (k) is the gluon propagator in the hard dense loop (HDL) approximation, (Gr(°)(p)) , is the inverse of the 

bare fermion propagator as is given in Eq. (g) above, and {G(p))~ b is the inverse of the full fermion propagator, defined 
over the true (color superconducting) vacuum. The latter, under the assumption of no wave function renormalization, 
is given by 



(G(p)Y 



1^7 



ii 



M7 



(4) 



The order parameter of the superconducting phase is represented by S, appearing in the off diagonal terms. 

Without loss of generality, in this paper we consider only the parity even order parameter. In the instanton approach, 
it is shown that such an order parameter determines the true vacuum of the color superconducting phase p,H . At 
sufficiently large density of quark matter when the instanton density vanishes (and the current quark masses are small 
compared to the scale of the condensate), the parity odd vacuum may have a relatively long life-time and, therefore, 
be of great interest too [p.6| . Obviously, the SD approach with the perturbative kernel, utilized here, is incapable of 
differentiating between the parity odd and parity even order parameters. The generalization of our analysis to the 
case when both parity odd and parity even contributions are non-zero is straightforward. 

Now, let us discuss the color-flavor structure of the order parameter. In Ref. 0, from symmetry arguments alone 
it is argued that the most general order parameter could contain two pieces: the color antitriplet, flavor antitriplet, 
^ab ~ fia^i ~ fia$b> an d tne c °l° r sextet, flavor sextet, S^ fc ~ 5 l a 8 b + S 3 a S l b . Therefore, the most general form of the 
solution for the order parameter reads £^ b ~ ni5 l a 8l + k 2 <5^ [jl4|. In this paper, we suggest an alternative (although, 
equivalent) representation of the color-flavor structure, 
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where 
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P 2 ^ \p\ 
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are the free quark (antiquark) on-shell projectors fgJlCJI, and P'"' are the new projectors in the color-flavor space, 

{v^Y ab ^ l -5 ab 5^- l -5i5l (7b) 

(^ (3) )1 = l 6ab6ij + \ & ° S * \*& (7c) 



These three projectors satisfy the condition of completeness: J2 n V^ = X where 2"* b = S^S 1 ^ . Also, all of them are 
symmetric under the simultaneous exchange of the color and flavor indices and, as a consequence of this property, 
(■p(«)) f = p(n) for n= 12) 3. 

As we mentioned above, the most general order parameter contains the antitriplct-antitriplet (3, 3) and the sextet- 
sextet (6,6) contributions. In addition to these, the three-parameter space swept by our projectors V^ n > contains 
an extra piece, the singlet-singlet S^ 6 ~ SabS 1 -* ■ In order to get rid of this singlet, it is sufficient to introduce the 



Then, as is easy to check, the two independent parameters, A* and A. 
,g 31 and the (6, 6) order parameter A, g 6 s 



restriction: A3 = — A2 

combinations of the (3, 3) order parameter A,„- *\ an d the (6, 6) order parameter A 



± 
2 ■ 



Af = 2 



A ± 

^(3,3) 



2A f 6 ,6) 



^ = A f3,3) 



A ± 

^(6,6)- 



are the linear 

(8) 
(9) 



(Note that our A ± and A 2 are related to K\ and K2 parameters of Ref. [H as follows: A t = 3k\ + K2 and A 2 = — K2O 
As we will see below, the advantage of working with the projectors V^ n ' is tremendous, while the price that one pays 
for such a convenience is minimal. 



By inverting the expression in Eq. (0), we obtain the propagator, 
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and €y = \p\ ± /u. In the derivation, we used the formulae (|Al|) - (A4) from Appendix |A|. 

Now, returning back to the SD equation (ph, we see that the sum over the adjoint indices in the right hand leads 
to the following result: 
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At this point, we recall that the color-flavor structure of the quark propagator is completely determined by that of 
the projectors in Eq. ([7|). Therefore, by using their definition, we derive the following relations: 
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By making use of the explicit form of the quark propagator in Eqs. (pOJ), (O) and the relations (fL3|), we derive the 
gap equation, 
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which, after taking into account the restriction A3 = —A* discussed earlier, splits into the following set of coupled 
equations: 
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So far, we did not say much about the form of the gluon propagator T>^„(q — p) used in the kernel of the SD 
equation. Apparently, the medium of high density quark matter modifies the bare gluon propagator in many different 
ways. The most important modification is given by the screening effects ]9|-^2[. By taking these into account in the 
HDL approximation, we arrive at the following expression for the gluon propagator |9| : 
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where k — q — p, M 2 — Nfa s fi 2 /it and the projectors 0\il are given by 
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These projectors define the three types of gluon modes: a~™ — O v ilv A v (magnetic), a^' = 2 V A" (electric) and 

ajl = OP v A v (longitudinal), respectively. The corresponding three terms in the right hand side of Eq. ( |l8| ) are the 
propagators of these three types of modes. 

Obviously, the gluon propagator is also modified by the Mcissner effect. However, as in the case of the Nf = 2 
QCD, this effect is negligible everywhere outside the narrow infrared region given by < p 2 < (Af 2 |Ao|) 2 ^ 3 |)] |l2| . 
By noting that the latter does not overlap with the most important dynamical region of color symmetry breaking, 
|Ao| 2 < p 2 < /i 2 , it is justified to neglect the Meissner effect while working at the leading order in the coupling 
constant. At the next to leading order, however, it will be crucial for the self-consistency of the solution to properly 
take this effect into account. 

Now, let us return back to the analysis of the gap equations in Eqs. ([16]) and (|l7J). We note that the main 
contribution in the right hand side of these equations comes from the terms proportional to A~ (we assume that /i is 
positive). Then, in our approximate analytical approach, it is justified to omit the A^~ terms altogether. As a result, 
the equations for A~ decouple, 
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Here we would like to emphasize that only Ar determine the values of the gaps in the spectrum of quasiparticles at 
the Fermi surface. In contrast, A+ just slightly modify the spectrum of anti-quasiparticles [see Eqs. (fLl|)], 

By following the same arguments as in Ref. |1J] , it can be shown that the dependence of the solutions for A~ (q) 
on the spatial momentum orientation is negligible. Then, the integral over the angular coordinates, parameterizing 
the relative orientation of q and p, can be calculated explicitly in the right hand side of the gap equations ( 20|) and 
(pl|). The results of these integrations for all types of gluon modes are presented separately in Eqs. ( |Aq ) - (A12) in 
Appendix |A|. 

Next, after neglecting the dependence of A^ on |<f| J9 H12J and approximating the kernels of the equations by their 
dominant contributions at the Fermi surface [see Eq. ( |A14| ) in Appendix 0] , we can also perform the integration over 
the absolute value of the spatial momentum. The result is presented in Eq. (A15) in Appendix [A]. By using it, we 
finally arrive at the following system of equations: 
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where v — y/8a s /9n and A = 16(2tt) 3 / 2 /j,/ (3a s ) 5 / 2 . This set of coupled equations is solved approximately in Ap- 
pendix |B[ Here we present the final result, 

A7(p 4 ) « 2A^(p 4 ) = -2Ai"(p 4 ) « 2A^ g) (p 4 ), (24) 

where 
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By recalling the definition in Eqs. (H) and rt7ft, we arrive at the following color-flavor structure of the order parameter: 

3 
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This is the dominant antitriplct-antitriplet portion of the solution for the order parameter in the color-flavor locked 
phase of QCD with three massless flavors. However, as one can see from the structure of the coupled equations 
jig) and (|17|), or even from their approximate version in Eqs. (^2|) and (p3|), the sextet-sextet portion of the order 
parameter, ATg 6) — (A^ — 2A^~)/6, should also be nonzero. 



Our analysis in Appendix |B| [see Eq. (B7) there] clearly demonstrates that the repulsion rather than attraction 
prevails in the channel with the sextet-sextet order parameter in the most important dynamical region A 2 , < p 2 < /i 2 . 
Therefore, a nonzero sextet-sextet contribution is generated only in the narrow infrared region < p 2 J$ A 2 , where the 
nonlinear character of the gap equations is essential. Outside this dynamical region, AT, e Jp) ~ []. Clearly, a more 
rigorous estimate of the magnitude of this sextet-sextet order parameter would be unreliable in the approximation we 
use. Indeed, in the infrared region of momenta, < p 2 < A 2 ,, all the subleading effects (the Meissner effect, the wave 
function renormalizations, the vertex corrections, the corrections from the Af terms), which we left aside, should be 
taken into account more carefully. 

Therefore, in further studies, it would be of great interest to analyze the gap equations up to the next to leading 
order. In addition to giving a proper estimate for the sextet-sextet contribution, the next to leading corrections could 
also produce a finite correction to the magnitude of the antitriplet-antitriplet order parameter. Indeed, the naive 
analysis of our gap equations, that ignores all the subleading corrections but the effect of the gaug e fixi ng term in the 



gluon propagator, would modify the value of the gap as follows: A,= g-, —> exp(3£/2)A f B ^ [see Eq. (A14)]. Apparently, 
the proper consideration of all the subleading effects would result in cancellation of any such a dependence on £. It is 
difficult, however, to show this explicitly, and it is not clear what kind of a correction this may produce for the value 
of the gap. 

In conclusion, in this paper we derived the coupled set of gap equations for the order parameters of the color-flavor 
locked phase of massless Nf = 3 dense QCD. Our analysis of these equations shows that the dominant contribution 
to the order parameter comes in the color antitriplet, flavor antitriplet channel. The other, the color sextet, flavor 
sextet contribution is small but non-zero. 

After this paper was completed, we learned about the paper by T. Schafer discussing somewhat similar issues [ pT[ . 
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APPENDIX A: FORMULAE USED IN DERIVATION OF GAP EQUATIONS 



In derivation of Eqs. (Jill) , we find it convenient to use the following representations for the kinetic terms of the 
spinor propagators: 
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where the projectors Ap are defined in Eq. (J5|) in the main text. 

In order to perform the angular integration in the right hand sides of the gap equations (16h and (J17|) , one first 
need to calculate the following two types of traces over the Dirac indices: 
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where £ = cos 6 is the cosine of the angle between three- vectors q and p, and irrelevant antisymmetric terms are 
denoted by ellipsis. 

By contracting these traces with the projectors of the magnetic, electric and longitudinal types of gluon modes, 
defined in Eqs. (|l9|), we arrive at 
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where q = |<f|, p = \p\, g 4 = -jgr and p 4 = -«Po- 

Now, by making use of the explicit expressions ( A7), we easily perform the angular integrations of all types appearing 
in the gap equations, 
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where u)f = (it/2)M 2 uj and u> = \qi — p 4 |. Similarly, 
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As a result, the angular average of the gluon propagator (multiplied by g 2 weight) in the vicinity of the Fermi 
surface is given by the following approximate expression: 
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where eZ = q — p and £ if the gauge fixing parameter. The three terms in the last expression are the leading 
contributions of the magnetic, electric and the longitudinal gluon modes, respectively. Because of the absence of 
the logarithmic factor in front of the gauge fixing parameter, the longitudinal gluon modes become relevant only in 
the next to leading order. Therefore, here we ignore it in what follows. (Notice, however, that the solution to the 
gap equation might still be very sensitive to the choice of the gauge in the approximation with no wave function 
renormalization taken into account.) Finally, the integral over the spatial momentum gives 
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q 2 dq 



o g 2 + ( e ,-) 2 + |A| 2 
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(2p) 3 



In 



(2/^ 



2Af 2 



Mix 2 



A 



3x/g 4 2 + |A| 2 



(A15) 



where A = {2^f / {y/2-wM 5 ) = 16(2^) 3 / 2 ^/(N f a s ) 5 / 2 



APPENDIX B: THE APPROXIMATE SOLUTION OF THE GAP EQUATIONS 



In order to solve the set of integral equations ( |22|) and (23), we further approximate the kernels and introduce finite 
infrared and ultraviolet cutoffs, 



. , . v 2 I f p dq A , , , A 
AiW - 7 / -A 2 g ln- 

2 V Ja q p 



d q a / m A 

— A 2 g ln- 

P q q 



(Bl) 



A . . v 2 ( f p dq A , ,. A [ A dq A , ^ A\ v 2 ( f p dq . . , , A f A dq 
A 2 (p)~- / ^A 1 (g)ln-+/ -*A ig lii- +— / -"±A 2 (g) In - + / -1 

16 \Ja q p Jp q q J » yj Ao q p J p q 



A 



A 2 (g)ln^j, (B2) 



where v — y/8a s /9ir and A = 16(27r) 3 /V/(3a s ) 5//2 . Also, for brevity of notation, we omit the superscript "-" at Aj 
in this appendix. The set of equations ( |Bl|) and (B2) is equivalent to the following two coupled differential equations: 



pA'{(p)+A' x {p) 
pA'>(p) + A' 2 (p) 



v 2 A 2 {p) 



0. 



A x (p) , v 2 A 2 {p) 



= 0, 



(B3) 

(B4) 



16 p 8 p 

along with the following infrared and ultraviolet boundary conditions: 

P A 'i»| p=Ao =0 (IR), A lj2 (A) = (UV). (B5) 

By introducing new variables, A/g^) and A( 66 ), as in Eqs. ([|) and @, the two differential equations decouple, 

H3,3)( 



P A (6,6)(p) + A '(6,6)(p) 



v 2 A fg . 3) (p) 



4 p 

V 2 A( 6 , 6 )(p) 

8 p 



= 0, 
= 0. 



(B6) 

(B7) 



Both the infrared and ultraviolet boundary conditions for A^agj and Ai ee \ are the same as those for Ai^(p), given in 



Eq. (B5). The negative sign in front of the v 2 term in Eq. (B7) indicates that the repulsion, rather than attraction, 
dominates in the pairing channel with the (6, 6) order parameters AjR^(p). This explains, in particular, why the only 
solution for A( 66 )(p), satisfying the boundary conditions in Eq. (B5), is the trivial one. The solution for A^g^p), 
on the other hand, reads 



A (3,3)b) = Aosinf-ln- j , A = Aexp(--J , 



(B8) 



where, in order to be consistent with the approximation used in derivation of Eqs. (Bl) and (B2), we defined the 
overall constant of the solution by requiring that A(g j g)(Ao) = Aq. 
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